RADIATIVE RETURN: A PROGRESS ON FSR TESTS 
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To improve the accuracy of the e+e~ — > 7r+7r~7 Radiative Return 
method, one has to control the theoretical uncertainty of the final-state 
photon emission. It is of particular importance at DAPHNE for the anal- 
ysis, where cuts are relaxed to cover the threshold region. By means of 
Monte Carlo generator PHOKHARA we compare several final-state radi- 
ation models and present results, relevant for a meson factory running at 
Vs = 1 GeV. 

PACS numbers: 14.40.Be, 13.40.Ks, 13.60.Le, 13.66.Bc, 13.66.Jn 

1. Introduction 

The Radiative Return Method [l,'2j 3] (RRM) allows an extraction of the 
hadronic cross section a^'^'^iCp) for hadronic invariant mass squared from 
the energy threshold up to the nominal energy of the experiment at the fixed 
beam energy e^e~ colliders. High-luminosity meson factories are especially 
suited for this purpose 0]. Interest in precise measurement of a^"''^[Q'^) 
is motivated, in part, by its relevance to the hadronic contribution to the 
muon anomalous magnetic moment a^^f^ [SJ [B] and the electromagnetic fine 
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structure constant [7j. However, the method can also be apphed to extract 
the meson form factors and other meson properties. 

The RRM uses da{e~^e~ — > hadrons+ photons) /dQ"^, a measured differ- 
ential cross section, for the extraction of da{e^e~ hadrons). For the the- 
oretical description, the perturbative QED diagrams at the leading order in 
QED coupling a (LO) and at the next to leading order (NLO) are considered 
and classified as initial-state radiation (ISR) or final-state radiation (FSR) 
ones. The kinematic cuts are applied to sufficiently suppress the FSR, when- 
ever possible, as the factorization d<7{s,Q )|/5^ = 

which allows for da^"''^{Q'^) extraction, holds for diagrams with ISR pho- 
tons only. The function TZ{s, Q^) is given by QED. The FSR part is model- 
dependent, thus dedicated numerical studies are needed for correct ISR- 
FSR separation. The Monte Carlo generator PHOKHARA was developed 
for these and related purposes: FSR at NLO has been included [8] for pion 
pair production and, in addition to scalar QED (sQED), some particular in- 
gredients (the (j) radiative decay) were implemented [9j. The FSR was also 
examined by other Monte Carlo programs, e.g., that with the Resonance 
Chiral Theory (R^T) motivated framework [lOJ and phenomenologically- 
oriented model which was also included into PHOKHARA 6.1 [12j . 

The reaction e~^e~ tt'^tt~j was explored by KLOE [13]: the cross sec- 
tion da^°'^/dQ'^ and pion form factor Ft,{Q'^) in the range 0.35 GeV^ < < 
0.95 GeV^ were extracted [H] from the on-peak {y/s = = 1.02 GeV) 
data sample by means of the RRM. However, the kinematic cuts, which 
were applied in order to suppress FSR, did not allow to measure at 
below 0.35 GeV^. 

One can measure the Ft^^Q"^) in the threshold region relaxing some of 
the cuts, but then one has to subtract the FSR contribution. In this sce- 
nario, one needs to control the description of the final-state emission process 
and detailed studies are needed to estimate the theoretical uncertainty. To 
simplify the analysis, it is better to perform the measurement off the (j) me- 
son peak, because in this case the contributions from the (p meson radiative 
decays are small and the FSR models can be controlled easier. 

The investigations presented here are of particular importance for the 
forthcoming KLOE RRM analysis of pion pair production. We focus on the 
off-(/)-peak measurement, at e"*"e~ center-of-mass energy ^/s = 1 GeV, for 
which KLOE collected 230 pb~^ of data [K]. Due to the interest in precision 
at small (i-e., below the p resonance), the Chiral Perturbation Theory 
(xPT) [16j can be helpful. The relevant theoretical aspects are sketched in 
Section [2j We use Monte Carlo generator PHOKHARA to compare several 
final-sate radiation models. The theoretical heritage of Virtual Compton 
Scattering (VCS) off the pion in xPT framework [171 [18] is used to estimate 
the role of higher order effects. 
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The numerical results for cross section and asymmetry are presented in 
Section [3l All the parameters of implemented models are fixed indepen- 
dently, thus one deals with model predictions. In Section [5] we present our 
conclusions. 



The transition 7* tt^tt^'j is described by the model-dependent FSR 
tensor M'^'ill- In all realistic models it contains the same Born-level contri- 
bution M^'^^^, which corresponds to a no-structure approximation for pion 
(Scalar QED or lowest-order xPT). Thus we consider -/Vf|^^^,^ as a model- 
independent part. 

The first correction accounts for the pion structure by means of the pion 
form factor. It replaces the Born-level amplitude by "Generalized Born" 
(GB) one [H], M^'"^, which is also called "sQED*VMD" |8]. Generalized 
Born FSR tensor reads 



where P and k are the virtual and real photon momenta, qi and 52 — pion 
momenta {Q = qi + q2)- 

Limit Ft^{P'^) —>■ 1 reproduces the M'^^^^ amplitude. Notice, that in 
e^e~ — > 7r"'"7r~7 at LO, = s, thus the form factor Ft^{P'^) 7^ 1 and 
its correction is never negligible. Also this part is well established both 
theoretically and experimentally. 

In the ISR amplitude, with 7* vr+vr" transition in the final state, one 
finds Ft^{Q'^) factor in the amplitude. Therefore, the tt+tt" invariant mass 
distribution is governed by the form factor shape. For consistency, one has 
to use the same expression for the pion form factor in the ISR and FSR am- 
plitudes. It is important to take the form factor tested experimentally and 
not to rely only on a particular model assumptions. This will be illustrated 
in the next Section. In order to understand the accuracy of approxi- 
mation, we study further corrections using the models of Refs. pT | [T7 1 118]. 



The first model, "VMD*xPT" , is based on 0{p'^) xPT SU{2) description 
of VCS 7*7r± ^ 77r± [17J and that in 5C/(3) case US]. The FSR tensor 
has the form M^^'^ = Mq'q + ^n'b- ^^^^ term is given by Eq. ([1]) 

and a straightforward improvement beyond xPT is supposed (denoted by 
prefix "VMD*"): the pion form factor is an external input (e.g., defined 
by parametrization of the measured Ft^). The second term, the Non-Born 



2. Theoretical issues of final-state radiation 



M^^^=-ie^F.(p2)(^ 



(fc + gi-g2)^gl ^ (k + q2-qiYq2 
qi-k q2-k 
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Fig. 1. Differential cross section for e+e^ — > 7r+7r^7. All the listed models, but 
the one with the 0(p*) xPT form factor, give very close predictions (overlapping 
curves in the plot) 



correction to FSR reads: M'^^ 
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(3) 



in SU (2) and SU (3) framework, correspondingly; see original papers pTtfTS] 
for the explicit form of the loop functions G-,^ and Gk- Numerical values of 
the low energy constants are F = 92.4 MeV, {Iq — I5) = 3.0 ± 0.3 and 
Fo = 87.7 MeV, as cited in [21], and {Ug + L\q) = (1.32 ± 0.14) x lO^^ ^t 
scale /i = Mp, as estimated in [22] . 

The second model [11], called the "main model" further in the text, 
can be considered as a parametrization of tt^tt^^ KLOE data, transformed 
to 7r"'"7r~7 via isospin symmetry [23]. It was implemented in FASTERD 
Monte Carlo generator [H] and in PHOKHARA 6.1 recently [12]. The FSR 
tensor contains Mq^^ given by Eq. ([T|) and the Non-Born corrections due to 
important vector-resonance and double-vector-resonance contributions. 



3. Numerical results 

We use Monte Carlo generator PHOKHARA to compare the model- 
dependent effects in e~^e~ — > tt~^tt~j cross section and asymmetry for the 
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Fig. 2. Non-Born corrections to da/dQ'^ (left) Note, that at NLO the interference 
is neglected. At > 0.5 GeV^ all the curves but "main model" overlap 



off-peak case at (?!)-meson factory, ^/s = 1 GeV. 

First of all, we stress that any simplification of the pion form factor 
can drastically affect the model results. Figm'e [1] shows that rigorous 0{p^) 
xPT form factor [16] gives completely wrong estimate for differential cross 
section even in the region of below the p meson peak. The theoretical 
explanation of the form factor role was given above. The form factor used 
in VMD*xPT and "main model" is the parametrization of available data 
given by Gounaris-Sakurai version of Ref. [20] . 

In Fig. [H one can see the very close cross section predictions, despite 
the fact, that the models have completely different Non-Born corrections. 
This is due to the fact that the GB contribution dominate for the given 
event selection. Taking the GB approximation, Eq. ([T]), as a reference, we 
plot {da[model] — da[GB])/da[GB]. To show the relative contribution of 
loop and "constant" terms in xPT we consider also the case of (/g — ^5) 
and (Lg -|- Liq) being artificially set to zero. Corresponding results are 
marked as "loop only" in the pictures. Figure [2] shows that the Non-Born 
corrections are at a few per cent level. From Fig. [3] one concludes that the 
FSR contribution to the cross section is significant in the whole range of 
Q^, especially at low Q^. 
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Fig. 3. The role of FSR in da/dQ^ (right). AU the hsted models give very close 
predictions (overlapping curves in the plot). Note, that at NLO the interference is 
neglected 



Pion forward-backward asymmetry (FBA) as a function of reads 

4 m+>90")-jV(g,+ <90°) 2^ 

^^^^^ ) = N{e^, > 90") + N{e^, < 90°) ^ 

in terms of numbers of events. Origin of the non-zero FBA is the interference 
of C-odd and C-even amplitudes, e.g., that of ISR and FSR at LO. Thus, 
FBA is sensitive to the relative phase, which may differ among the models 
even if they predict the same cross section. Notice, that the experimen- 
tal data on asymmetry and cross section are to large extent independent. 
Therefore the FBA is a good test for models. Aspects of using the FBA in 
e+e~ —!' 7r+7r~7 were discussed in O El [9l [TOl [19]. 

Figure m shows that FBA is sizable and relatively easy measurable. From 
Fig. Owe conclude that the Non-Born corrections to the FBA are of few per 
cent order and will not have a big influence on the theoretical uncertainty. 

It has to be stressed that if the xPT corrections were not accounted for 
in the formulae used to measure the li'^Trl, they are partly accounted for in 
the experimental parameters of F.,^ and other model parameters. In other 
words, one model should be used in all experimental analyzes and adding 
ad hoc additional corrections is not appropriate. 

4. Conclusions 

Using PHOKHARA, we studied the corrections given by [17\ I18j. 
and by a phenomenological model including miscellaneous hadronic reso- 
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Fig. 4. Forward-backward asymmetry for pion. All the listed models give very close 
predictions (overlapping curves in the plot) 



nance effects [H]. Corrections due to ai resonance [19] are to be considered 
elsewhere. The role of the pion form factor is seen to be very important. 
Final-state radiation is significant in the whole range of Q^, especially at 
low Q^. We have found the NLO corrections to be non-negligible, even if 
the Generalized Born contribution is dominant. Non-Born corrections are 
of order of few per cent. They diff'er among the models, but it will be dif- 
ficult to distinguish them with the present KLOE ofi'-peak statistics. The 
results presented here show that one should include the xPT corrections in 
the analysis when the accuracy of the experiment reaches a per cent level. 
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